Microgeographic evolution of metabolic physiology in a salamander metapopulation.

The Metabolic Theory of Ecology explains ecological variation spanning taxonomic organization, space, and time based on universal physiological relationships. The theory depends on two core parameters: the normalization constant, a mass-independent measure of metabolic rate expected to be invariant among similar species, and the scaling coefficient, a measure of metabolic change with body mass commonly assumed to follow the universal 3/4 scaling law. However, emerging evidence for adaptive microevolution of metabolic rates led us to hypothesize that metabolic rate might exhibit evolved variation among populations on microgeographic scales. To evaluate our hypothesis, we explored evidence for evolved variation in the scaling coefficient and normalization constant within a spotted salamander (Ambystoma maculatum) metapopulation in Connecticut, USA. We measured standard metabolic rate in common-garden raised spotted salamanders from 22 different populations and tested for the effects of six ecological variables suspected in advance to select for divergent physiology. We found that metabolic rate rose with body mass with a log-log slope of 0.97 that was statistically different from the expected 3/4 scaling law. Although we found no evidence for interpopulation variation in the scaling coefficient, we found evidence for interpopulation variation in the normalization constants among populations. Metabolic variation was best explained by differences in population density among ponds. Our results provide mixed support for Metabolic Theory of Ecology assumptions about parameter invariance and illustrate how fundamental physiological processes such as metabolic rate can evolve across microgeographic spatial scales.

Balancing selection and drift in a polymorphic salamander metapopulation.

Understanding how genetic variation is maintained in a metapopulation is a longstanding problem in evolutionary biology. Historical resurveys of polymorphisms have offered efficient insights about evolutionary mechanisms, but are often conducted on single, large populations, neglecting the more comprehensive view afforded by considering all populations in a metapopulation. Here, we resurveyed a metapopulation of spotted salamanders (Ambystoma maculatum) to understand the evolutionary drivers of frequency variation in an egg mass colour polymorphism. We found that this metapopulation was demographically, phenotypically and environmentally stable over the last three decades. However, further analysis revealed evidence for two modes of evolution in this metapopulation-genetic drift and balancing selection. Although we cannot identify the balancing mechanism from these data, our findings present a clear view of contemporary evolution in colour morph frequency and demonstrate the importance of metapopulation-scale studies for capturing a broad range of evolutionary dynamics.

Direct effects influence larval salamander size and density more than indirect effects.

Direct and indirect effects both influence population and community dynamics. The relative strengths of these pathways are often compared using experimental approaches, but their evaluation in situ has been less frequent. We examined how individual and aggregate impacts of direct and indirect effects of species densities, proxies for competition and predation pressure, and habitat variables influenced patterns of larval density and body size of ringed (Ambystoma annulatum) and spotted salamanders (A. maculatum). We surveyed > 150 ponds in Missouri, USA, from 2012 to 2014 to measure the density and body size of each focal species, the density of co-occurring pond food web members, and select habitat features. We used structural equation modeling to quantify the relative importance of direct and indirect pathways on both body size and larval density. Overall, both responses were explained through a combination of direct and indirect effects. However, the magnitudes of direct effects were often greater than indirect effects. Some of the direct and indirect relationships with larval salamander size and density were also consistent with results from experimental studies. Finally, total direct and indirect effects were often weaker due to habitat and density variables negating each other's impacts. Overall, our study shows that direct effects were equivalent to, or more important than, indirect effects. We also demonstrate that the effects stemming from individual relationships can sum to produce net patterns that are negligible in magnitude. Further work on direct and indirect effects with observational data are needed to examine their magnitudes in natural communities.

The influence of breeding phenology on the genetic structure of four pond-breeding salamanders.

Understanding metapopulation dynamics requires knowledge about local population dynamics and movement in both space and time. Most genetic metapopulation studies use one or two study species across the same landscape to infer population dynamics; however, using multiple co-occurring species allows for testing of hypotheses related to different life history strategies. We used genetic data to study dispersal, as measured by gene flow, in three ambystomatid salamanders (Ambystoma annulatum, A. maculatum, and A. opacum) and the Central Newt (Notophthalmus viridescens louisianensis) on the same landscape in Missouri, USA. While all four salamander species are forest dependent organisms that require fishless ponds to reproduce, they differ in breeding phenology and spatial distribution on the landscape. We use these differences in life history and distribution to address the following questions: (1) Are there species-level differences in the observed patterns of genetic diversity and genetic structure? and (2) Is dispersal influenced by landscape resistance? We detected two genetic clusters in A. annulatum and A. opacum on our landscape; both species breed in the fall and larvae overwinter in ponds. In contrast, no structure was evident in A. maculatum and N. v. louisianensis, species that breed during the spring. Tests for isolation by distance were significant for the three ambystomatids but not for N. v. louisianensis. Landscape resistance also contributed to genetic differentiation for all four species. Our results suggest species-level differences in dispersal ability and breeding phenology are driving observed patterns of genetic differentiation. From an evolutionary standpoint, the observed differences in dispersal distances and genetic structure between fall breeding and spring breeding species may be a result of the trade-off between larval period length and size at metamorphosis which in turn may influence the long-term viability of the metapopulation. Thus, it is important to consider life history differences among closely related and ecologically similar species when making management decisions.

Life history differences influence the impacts of drought on two pond-breeding salamanders.

Drought is a strong density-independent environmental filter that contributes to population regulation and other ecological processes. Not all species respond similarly to drought, and the overall impacts can vary depending on life histories. Such differences can necessitate management strategies that incorporate information on individual species to maximize conservation success. We report the effects of a short-term drought on occupancy and reproductive success of two pond-breeding salamanders that differ in breeding phenology (fall vs. spring breeder) across an active military base landscape in Missouri, USA: We surveyed ~200 ponds for the presence of eggs, larvae, and metamorphs from 2011 to 2013. This period coincided with before, during, and after a severe drought that occurred in 2012. The two species showed contrasting responses to drought, where high reproductive failure (34% of ponds) was observed for the spring breeder during a single drought year. Alternatively, the fall breeder only showed a cumulative 8% failure over two years. The number of breeding ponds available for use in the fall decreased during the drought due to pond drying and/or a lack of re-filling. Estimates of occupancy probability declined for the fall-breeding salamander between 2012 and 2013, whereas occupancy probability estimates of the spring breeder increased post-drought. The presence of fish, hydroperiod, the amount of forest cover surrounding ponds, and canopy cover were all found to affect estimates of occupancy probabilities of each species. Pond clustering (distance to nearest pond and the number of ponds within close proximity), hydroperiod, forest cover, and canopy cover influenced both estimates of colonization and extinction probabilities. Our results show life history variation can be important in determining the relative susceptibility of a species to drought conditions, and that sympatric species experiencing the same environmental conditions can respond differently. Consideration of the spatial network and configuration of habitat patches that act as refuges under extreme environmental conditions will improve conservation efforts, such as the placement of permanent ponds for aquatic organisms. A better awareness of species-specific tolerances to environmental filters such as drought can lead to improved management recommendations to conserve and promote habitat for a greater diversity of species across landscapes of spatially connected populations.

Intermediate pond sizes contain the highest density, richness, and diversity of pond-breeding amphibians.

We present data on amphibian density, species richness, and diversity from a 7140-ha area consisting of 200 ponds in the Midwestern U.S. that represents most of the possible lentic aquatic breeding habitats common in this region. Our study includes all possible breeding sites with natural and anthropogenic disturbance processes that can be missing from studies where sampling intensity is low, sample area is small, or partial disturbance gradients are sampled. We tested whether pond area was a significant predictor of density, species richness, and diversity of amphibians and if values peaked at intermediate pond areas. We found that in all cases a quadratic model fit our data significantly better than a linear model. Because small ponds have a high probability of pond drying and large ponds have a high probability of fish colonization and accumulation of invertebrate predators, drying and predation may be two mechanisms driving the peak of density and diversity towards intermediate values of pond size. We also found that not all intermediate sized ponds produced many larvae; in fact, some had low amphibian density, richness, and diversity. Further analyses of the subset of ponds represented in the peak of the area distribution showed that fish, hydroperiod, invertebrate density, and canopy are additional factors that drive density, richness and diversity of ponds up or down, when extremely small or large ponds are eliminated. Our results indicate that fishless ponds at intermediate sizes are more diverse, produce more larvae, and have greater potential to recruit juveniles into adult populations of most species sampled. Further, hylid and chorus frogs are found predictably more often in ephemeral ponds whereas bullfrogs, green frogs, and cricket frogs are found most often in permanent ponds with fish. Our data increase understanding of what factors structure and maintain amphibian diversity across large landscapes.

Habitat traits and species interactions differentially affect abundance and body size in pond-breeding amphibians.

In recent studies, habitat traits have emerged as stronger predictors of species occupancy, abundance, richness and diversity than competition. However, in many cases, it remains unclear whether habitat also mediates processes more subtle than competitive exclusion, such as growth, or whether intra- and interspecific interactions among individuals of different species may be better predictors of size. To test whether habitat traits are a stronger predictor of abundance and body size than intra- and interspecific interactions, we measured the density and body size of three species of larval salamanders in 192 ponds across a landscape. We found that the density of larvae was best predicted by models that included habitat features, while models incorporating interactions among individuals of different species best explained the body size of larvae. Additionally, we found a positive relationship between focal species density and congener density, while focal species body size was negatively related to congener density. We posit that salamander larvae may not experience competitive exclusion and thus reduced densities, but instead compensate for increased competition behaviourally (e.g. reduced foraging), resulting in decreased growth. The discrepancy between larval density and body size, a strong predictor of fitness in this system, also highlights a potential shortcoming in using density or abundance as a metric of habitat quality or population health.